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SUMMARY 


This report describes the formulate cn and Implementation of 
navigation systems used for research Investigations In the NASA 
Ames V/STOLAND avionics system. The navigation systems described 
provide position and velocity in a cartestlan reference frame 
aligned with the runway. They use filtering techniques to com- 
bine the raw position data from navaids (e.g., TACAN, MLS) with 
data from onboard inertial sensors. The inertial data can be 
from either a low quality strapped down system or a precision 
platform inertial navigation system. The filtering techniques 
described use both complementary and Kalman filters. 

This report describes the software for the navigation systems. 
Later reports will provide the evaluation and analysis of the 
performance of these systems. 
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INTRODUCTION 


This report describes the formulation and implementation of navigation 
systems used for research investigations in the NASA Ames V/STOLAND avionics 
system. The V/STOLAND system was developed by Sperry Flight Systems and is 
an all-digital avionics system implemented in a UHIH helicopter. The system 
can operate using manual, stability augmentation or fully automatic modes 
in controlling the helicopter. The overall system capabilities permit re- 
search investigations for the development of all-weather capability for 
VTOL vehicles operating from small heliports. 

The VTOL aircraft approach, landing and departure paths can differ 
significantly from those of conventional aircraft at an air terminal. For 
example , the VTOL descent path can use a helix to avoid airspace used by 
the conventional aircraft. Following descent, a deceleration maneuver to 
a hover condition over the small landing pad can be made. Following hover, 
a letdown maneuver to the pad can be executed. In order to perform such 
maneuvering and landing under all-weather conditions the navigation system 
must produce accurate estimates of the aircraft position and velocity. 

The V/STOLAND system is currently implemented with the following 
navigation equipment. 

1. Microwave landing system (MLS) receiver providing range, 
azimuth and elevation measurements. 

2. TACAN receiver providing range and bearing measurements. 

3. VOR/DME receiver providing range and bearing measurements. 

4. Barometric altitude. 

5. JTEC Instrument for measuring true airspeed. 

6. R.i«dlo altimeter. 

7. Vertical gyro for aircraft pitch and roll measurements. 
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8. Directional gyro for aircraft heading measurement. 

9. Triad of body mounted accelerometers. 

10. A platform type inertial navigation system (LTN 51). 

All measurements are available in the dual Sperry 181 9B computers of the 
V/STOLAND system. 

Flipht tests of the system are conducted at the NASA/Ames flight test 
facility at Crows Landing, California. This facility is equipped with the 
Microwave Landing System (MLS) as well as radar tracking equipment, recording 
equipment and other features which provide the capability of evaluation of 
navigation systems which use the ’ii S as the primary navigation reference. 

The simulated heliport at the facility is a small area on the runway located 
approximately 2100 meters fror, the MLS range and azimuth station and about 
1000 meters from the MLS elevation station. 

The navigation systems described in this report all provide position and 
velocity in a cartesian reference frame with the "1" axis (x axis) aligned 
with the runway (at Crows Landing) and the "3" axis (z axis) aligned with 
the local vertical. TACAN data and barometric altitude provide position 
information in the terminal area prior to acquiring MLS. MLS range and 
azimuth are used through touchdown for the x>-y position information. MLS 
elevation data are used for vertical information when within the coverage 
area to an altitude of about 152 m. Thereafter radio altimeter data are 
used for the altitude information. 

The navigation systems use filtering techniques implemented in the 
airborne computer software to combine the raw position data with data 
from inertial reference Information from either (a) the body mounted 
accelerometers and the vertical and directional gyros (strapped- down IMlO 
or (b) Inertial navigation system (INS). The navigation systems described 
herein use the following filter types and inertial information. 
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1. A complementary filter using the body-mounted Inertial Information. 

2. A complementary filter using the LTN-51 inertial information. 

3. A Kalman Filter using the LTN-51 inertial information. 

4. A Kalman Filter using the body-mounted inertial information. 

These navigation systems are being flight tested to demonstrate their 
accuracy while varying the software's complexity and the accuracy of the 
onboard inertial-sensing equipment. 


The objective of this report is to describe the software for the 
navigation systems. Later reports will provide the evaluation and analysis 
of the performance of these systems. 
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NOTATION AND DEFINITIONS 


Notation 


The notation of over a symbol has the customary meaning of 
differentiation with respect to time. The (hat) mark over a symbol 
means the "estimated" or "computed" value of the symbolized quantity. 

The letter "d" before a quantity indicates an error or small variation of 
that quantity. For example, if X is the true value of position, it may 
be written as the sum of the estimated position and the position error, or 

X - X + dx . 


time. 


The notstion t, , t, , 
k’ k+1 

The time point 


, etc., are used to denote discrete points in 
occurs At seconds after tj^, or 


The time increment At denotes the primary cycle time of the implemented 
digital filter. 


Roman Symbols 


A, A 


b ,b .b 

ax ay az 

AAA 

b ,b ,b 
ax ay* az 




c ,c 

X V 

dx 


dx 

dx,dy,dz 

dx.dy.dz 

dz 


e 

F 

X 

g 


H.H 

m 


H 


e 


ma 


mr 


h 

\b 

^B’^R 


ma 


‘RW 


- discrete form of the matrix F . 

X 

-* acceleration measurement biases. 

- estimates of acceleration measurement biases. 

• actual and estimated bias error li< the TACAN range 
measurement . 

- actual and estimated bias error in the TACAN bearing 
measurement . 

- position and velocity smoothing vectors. 

A 

- the continuous error state vector of the estimate X for 
the level Kalman filter. 

- filter estimate of the error state vector dx. 

- errors in x,y,z components of position 

- errors lx x,y,z components of velocity 

A 

- the continuous error state vector of the estimate z for 
the vertical Kalman filter 

- base of natural logarithms, e « 2.7182... 

- n X n system dynamics matrix. 

- acceleration due to gravity 

- exter al measurement distribution (sensitivity) matrix. 

- sensitivity vector of altitude derived from the MLS 
elevation measurement to entimated state. 

- sensitivity vector of MLS azimuth to estimated state. 

- sensitivity vector of MLS range to estimated state. 

- estimated altitude in vertical complementary filter. 

- bias in the barometric altitude measurement 

- baro altitude and radio altitude. 

- altitude above ground measured by the radio altimeter 
or MLS elevation. 

- altitude of aircraft above MLS range/azimuth station. 

- raw altitude from complementary filter's measurement 
selection logic. 

• runway altitude with respect to sea level. 
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I 

Kj.Kj.K 


L ,L ,L 
X y z 

L* ,L* ,L* 

X y z 


Q 

Qe 





R 

e 



R 


T 


R 

V 




r 

c 

s 

t 


X,X 



- identity natrix. 

- feedback gains in the vertical complementary filter. 

- limit values on magnitudes of position and velocity 
smoothing vectors. 

•• magnetic north. 

- assumed variance of the random error q. 

- assumed variance of MLS (elevation derived) altitude 
measurement. 

- assumed variance of MLS azimuth measurement noise. 

- asstuned variance ot MLS range measurement noise. 

- random noise error in the external measurement. 

- random noise in the MLS (elevation derived) altitude 
measurement 

random noise error in MLS azimuth measurement, 
random noise in MLS range measurement. 

- estimated slant range to MLS elevation antenna. 

- estimated ground range. 

TACAN measured slant range. 

- VOR/DME measured slant range 

- MLS elevation station slant range ard its component 
in the x-y plane. 

- TACAN or MLS measured ground range. 

- Laplace transform variable. 

- t ime . 

- actual and estimated values of the aircraft's state vector. 

- x-component of aircraft velocity relative to the airmass. 

- x-component of acceleration from IMU. 

- navaid-derived position from complementary prefilter. 

- x-component of wind velocity. 
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< >> 


r 


x.y.* 


*b 


z 

e 



y 


a 


position of aircraft in a Cartesian reference frame 
with X along the runway, y in the horizontal plane 
normal to the runway, and z normal to the horizontal 
plane and positive pointing downward. 

filter estimates of x,y,z. 

vector of aircraft acceleration in body coordinates. 

location components of the MLS elevation antenna with 
respect to the runway reference frame. 

vector of aircraft acceleration in inertial platform 
(LTN-51) coordinates. 

aircraft position coordinates with respect to the MLS 
elevation antenna. 

horizontal true airspeed components in the runway 
reference frame. 


*RA’^RA’*RA 

*R’^R’*R 


X and y (runway coordinates) accelerations measured by IMU. 

location components of the MLS range and azimuth station 
with respect to the runway reference frame. 

aircraft position coordinates with respect to the run- 
way reference frame. 

raw acceleration in the runway reference frame as com- 
puted by V/STOLAND system software. 

estimated position vector. 




X ,y ,z 

V v V 


y 

e 



y 


location components of the TACAM ground station with j 

respect to the runway reference frame. ] 

location components of the VOR/DME ground station with | 

respect to the rnnway reference frame. 

computed value of the external measurement. j 

I 

actual and estimated MLS elevation-derived altitude measurements. i 

i 

external state measurement of aircraft. I 


m 


y ,y 

ma ma 

actual 

and 

A 

Y ,y 

mr mr 

- actual 

and 


estimated MLS azimuth measurements, 
estimated MLS range measurements. 


i 
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i 



major time (cycle) update of the Kalman filter, 
acceleration integration period of the filter. 

MLS elevation angle measurement, 
summation 

standard deviation (std) of acceleration colored 
noise (x,y,z channels). 

std of MLS elevation random noise. 

std of barometric altitude colored noise (bias). 

std of MLS azimuth random noise. 

std of MLS range random noise. 

std of radio altimeter random noise. 

std of TACAN bearing random noise. 

std of TAC\N '’:.iige random noise. 

std of TACaN range colored noise (bias). 

std of velocity noise. 

std of air data velocity noise. 

std of wind colored noise (bias) components. 

std of TACAN bearing colored noise (bias), 
std of initial heading measurement. 

time constant of an exponential function. 

time constant for acceleration measurement colored 
noise (x and y channels). 

time constant for TACAN range measurement colored noise. 



time constants used to compute smoothing vectors 
c and c . 

X V 

time constant for TACAN bearing measurement colored 
noise. 
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*u<Vl"k> 


Ip, 6,4) 


'^r’'^RW 


‘"lx’‘"2x’^3x 


- state transition matrix from time point t^^ to 
time point 

- forcing function sensitivity matrix affecting 

state ac due to . 

- aircraft attitude angles (roll, pitch, heading) 
measured by vertical and directional gyros. 

- bearing of runway with respect to magnetic north. 

- TACAN measured bearing from magnetic north. 

- VOR measured bearing from magnetic north. 

feedback gains in the complementary x filter. 


Abbreviations and Acronyms 


a/c 

- 

aircraft. 

ADI 

- 

attitude director indicator. 

AGL 

- 

above ground level. 

AUTO 

- 

automatic navigation mode of the V/STOLAND system. 

brg 

- 

bearing. 

rss 

- 

control stick steering mode of the V/STOLAND system. 

CTOL 

- 

conventional take-off and landing aircraft. 

DDAS 

- 

digital data acquisition system. 

HSI 

- 

horizontal situation indicator. 

Hz 

- 

Hertz, a unit meaning one cyclt ;sr second. 

IMU 

- 

inertial measuring unit. 

INS 

- 

inertial navigation system. 

I/O 

- 

input /output. 

MFD 

- 

multifunction display, a CRT on the research pilot's 
panel in the UH-lH. 

MLS 

- 

microwave landing system, provides range, azimuth and 
elevation. 

MODILS 

- 

modified instrument landing system, a precursor of MLS. 

MSL 

- 

mean sea level. 

MSP 

- 

mode select panel, located on the research pilot's 
panel in the UH-IH. 
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navald 

- 

an aid to aircraft navigation. 

RMDU 

- 

multiplexer and demultiplexer unit. 

rms 

- 

root mean square. 

mg 

- 

range . 

std 

- 

standard deviation. 

STOL 

- 

short take-off and landing aircraft. 

STOLAND 


automatic guidance, navigation and display system 
for short take-off and landing aircraft built by 
Sperry Flight Systems. 

TACAN 

- 

ultra-high frequency tactical air navigation aid, 
provides range and bearing information. 

VOR/DME 


very high frequency omni-range station with distance 
measuring equipment, which together provide range 
and bearing Information. 

V/STOLAND 


automatic guidance, navigation and display system 
for vertical and/or short take-off and landing air- 
craft built by Sperry Flight Systems. 

VTOL 

- 

vertlca] take-off and landing aircraft. 

WPT 

- 

navigation waypoint. 


TEST SYSTEMS 


This section gives a description of the V/STOLAND avionics system, 
the Crows Landing test facilities, the test flight paths used at Crows 
Landing and a sample of the navaid data for a typical approach. 


V/STOLAND System 

V/STOLAND is an integrated digital avionics system that provides 
navigation, guidance, control and displays for the UH-IH helicopter. For 
a good overall description of the V/STOLAND system see Ref. 1. The summary 
given here was obtained from Ref. 1. V/STOLAND is a flexible system, 
allowing evaluation of the aircraft's performance in various configurations 
of automatic control, display and navigation. The system is self-monltor- 
ing, with provisions for automatic disengagement when failures are de- 
tected by the system monitors. 

The V/STOLAND system provides the capability to fly conventional modes 
such as airspeed select and hold, altitude select and hold, flight path 
angle select and hold, and TACAN or VOR radial guidance modes. The Way- 
point (WPT) guidance mode also provides radial guidance to an arbitrary 
waypoint selected by the pilot. Conventional approaches are possible using 
selectable MLS glideslopes. With MLS navigation data, helical descent 
may be made to the touchdown point using a 334m radius 3-turn helix 
(HELIX mode) or using a 530m radius 2-tum helix (OFFSET HELIX mode). 

The MLS azimuth and rr.rge data are used for computing the deviation in- 
formation required for guidance. The system also provides for capturing 
and tracking a 3-dimensional reference flight path. In the basic navigation 
system, the TACAN, VOR/DME or MLS navigation sources may be selected 
manually or automatically fay priority logic which selects the most accurate 
and valid navaid. 
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V/STOLAND may be operated in three basic contrrl configurations, with 
or without the Flight Director: 

1 . Manual 

2. CSS (Control-Stick Steering) 

3. AUTO (Autopilot) 

In the manual configuration, the pilot controls the helicopter manually by 
the sticks and the pedals. No servos are engaged. If CSS is engaged, the 
Research Stick (left side) operates in a fly-by-wire mode, providing con- 
trol of the helicopter through the servos. If AUTO is engaged, the guid- 
ance and control of the system is fully automatic. 

When the Flight Director or AUTO is not engaged, conventional angular 
VOR/TACAN radial deviation Information is provided on the horizontal situa- 
tion indicator (HSI) , permitting manual capture and tracking of a selected 
radial. When the Flight Director is engaged, pitch, roll and collective 
commands are displayed on the attitude director indicator (ADI) . Linear 
rather than angular deviation from the course is always displayed on the 
HSI for any of the various guidance modes that may be selected on the Mode 
Select Panel (MSP) . 

The flexibility of the V/STOLAND system is significantly increased by 
providing for research modes which function through the Research Computer. 
Navigation, guidance, control and display modes may be programmed by the 
researcher in the Research Computer. The research modes can be exercised 
in any of the basic control configurations of the system. 

A block diagram of the V/STOLAND system is shown in Figure 1. The 
1819B is a general-purpose digital computer with a 18K memory of 18-bit 
words and capable of real-time operations in an airborne environment. 

The I/O is fully buffered and parallel and uses a party-line transmission 
system. The data adapter provides the required interface between the basic 
computer and the rest of the system. It performs all the analog-to-digital 
and digital-to-analog conversions and all digital-to-digital data transfers 
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The displays such as the ADI, HSI, and the MFD (multifunction display) 
provide the inertial, navigational and guidance information. The ADI 
includes the flight director command bars in addition to basic attitude 
data. The HSI provides navigation and guidance data. The MFD displays 
the horizontal situation of the aircraft and pertinent background data such 
as geographical features, navaid descriptors and the reference flight path. 


The servo complement consists of four series and four parallel servos. 
The series servos are electrohydraullc and have limited authority. Their 
movements are not reflected in the control sticks. The parallel servos 
are electromechanical rate servos that have nearly full authority, but 
limited rate capability, and result in movement of the sticks and pedals. 

The V/STOLAND system also interfaces with the Digital Data Acquisition 
System (DDAS) which consists of a Multiplexer and Demultiplexer Unit (RMDU), 
a tape recorder and a telemetry transmitter. The sensor input data as 
well as the computed data are transmitted to the DDAS for recording on 
tape. The taped data can be converted on the ground into strip-chart 
recordings for flight analysis. 

The emphasis of this report is on using the V/STOLAND system as a 
test bed for research on VTOL navigation systems which include the com- 
plementary filter navigation system located in the Basic Computer and the 
Kalman filter navigation system located in the Research Computer. 
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Crows Landing Test Facility 


Flight tests of the V/STOLAND avionics system are performed at the U.S. 
Navy’s Crows Landing Auxiliary Landing Field near Patterson, Cal. Figure 2 
shows a plan view of the field and locations of the TACAN station, the MLS 
range/azimuth antenna, and the MLS elevation antenna. One navaid not shown 
but occasionally used is the Modesto VOR/DME, which is 28 km NNE. 

All landing approaches are made to Runway 35. The runway coordinate 
system is defined such that the x-axis is on the runway centerline and 
positive northward, the z-axis is vertical and positive down, and the 
y-axis completes a right-handed Cartesian system. The origin is on the 
runway centerline and such that the MLS elevation antenna is in the y-z 
plane. The target touchdown point is at x ■ -914m, which is between 
Runway 35 's threshhold and the intersection of Runways 35 and 30. » 

In order to evaluate the performance of the V/STOLAND guidance and 
navigation, key variables are telemetered from the aircraft and recorded 
on the ground. Also, twin radars track the aircraft, and its position and 
velocity are also recorded. 


Test Flight Path 

The reference test flight path and approaches are shown in Figure 3. 
Three landing approaches have been programmed into the system. These are: 
(1) a straight-in approach of constant descending flight-path angle to a 
final 2.5° glideslope, (2) a three-turn helical descent with a radius of 
354m, and (3) a two-tum helical descent with a radius of 530m. The three- 
turn helix is not shown in Figure 3, but it begins and ends at the same 
points as the two-tum helix. Because of its greater radius, the two-turn 
helix requires a lower bank angle (about 10°) to fly and is more confortable 
for the pilots. Its main advantage, however, is that the bank angle 
excursions in windy conditions are considerably reduced. The glideslope 
for both helices is 6.11 deg. 
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After exit from either helix, the glideslope remains unchanged until 
the path intercepts a final glideslope of 2.5°(see Fig. 4). Finally, a 
constant-altitude segment is captured and flovm until hover is achieved 
over the touchdown point, and a short vertical "letdown" phase completes 
the landing. 


Sources of Navigation Information 


General characteristics of the sources of navigation information are 
discussed in this section. For this purpose telemetry data from a landing 
approach (the sixth of UH-IH Flight 9101) which flew the two-tum helix 
have been selected as typical. Figure 5 shows plots of the data from a 
time shortly before entry to the helix until just before letdown. 

The geometry of navigation with respect to either the TACAN station or 
VOR/DME station is shown in Figure 6. Slant range and magnetic bearing 
from the station are the measured quantities. Figure 7 depicts the geo- 
metry of MLS navigation, where the measured quantities are slant range, 
azimuth referenced to the runway centerline, and elevation angle. When 
available, MLS is used rather than TACAN, and TACAN is used rather than 
VOR/DME. 

As can be seen in Fig. 5a, the TACAN range measurement is fairly 
coarse; one standard deviation of the apparent noise on the signal is 
about 15 m. The recorded TACAN bearing signal (the second plot of Fig. 5a) 
has several spikes and discorilinuities due to a failure in the receiver 
which shifted the measured bearing by + 40°. The sinusoidal behavior of 
these signals (and many others in Fig. 5) between 90 and 290 seconds re- 
flects the helical flight path followed in the approach. The slant range 
in the second turn of the helix is slightly less than at a corresponding 
point in the first turn mostly because the altitude is less. The TACAN 
bearing signal improves markedly in the second turn of the helix because 
the elevation of the aircraft above the TACAN transmitter is substantially 
less than at previous times shown in the plot. This is entirely expected 
because the TACAN system is not designed for accuracy at high elevation 
angles. 
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It is clear from these plots that horizontal navigation for an approach 
and landing using TACAN data alone would be difficult. The microwave land- 
ing system (MLS) is much more accurate since it was designed specifically 
for this task. Therefore, it is logical that TACAN data are ignored if MLS 
data are available. 

Fig. 5A shows that the MLS range data are smooth but the azimuth signal 
has discemable aberrations, particularly between 90 and 140 seconds and 
between 200 and 220 seconds. These aberrations appear in virtually the 
same spot in nearly all helical approaches and are suspected to be due 
partly to the high elevation angle and partly to the signal's reception 
being switched from the fore to the aft MLS antenna on the UH-IH near 
these points of the helix. 

The MLS elevation data are obviously invalid between 70 and 130 seconds. 
Although the signal is smooth and appears valid after 130 seconds for this 
approach, it was not always reliable above 10° elevation. Therefore, the 
signal is ignored until the elevation angle, as computed from baro altitude 
and the estimated x and y positions, is below 10°. 

Figure 5b shows plots of the roll angle and pitch angle outputs of 
the vertical gyro, the heading angle output of the directional gyro, the 
true airspeed output of the JTEC sensor, and the outputs of the barometric 
and radio altimeters. As can be seen, the average roll angle in the helix 
is about 11° with variations to either side of about 3°. A small part of 
this variation is the result of flying a circular ground track in a steady 
wind; the bank angle is steeper when the aircraft is flying downwind and 
shallower when flying upwind. Because the reported wind was from 260 at 
2.6 m/s, the roll angle should be slightly steeper at a heading of 80° 
and shallower at 260° (-100°). However, it is clear from the plot that 
other significant aircraft maneuvering is taking place as well. 

The aircraft maintains about 2° pitch-up attitude throughout the helix 
with variations of 1° to 2° to either side. There is a steady increase in 
pitch attitude as the aircraft exits the helix, decelerates and flares. 
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True airspeed is not used by the Kalman filters except during initializa- 
tion, but it is used by the complementary filters for dead reckoning and 
wind estimation. The approach shown in Fig. 5b was flown at about 32 m/s 
(62 knots true), and this speed was maintained until exit from the helix. 

The airspeed measurement here shows greater variation than the 0.6 to 
1.5 m/sec peak-to-peak variation usually seen in this measurement. This 
may have been the result of some small turbulence aloft. 

The plot of the barometric altimeter reading shows that the approach 
begins at 805m (2640 ft) above mean sea level (MSL), which means the 
pattern altitude is 762m (2500 ft) above ground level (AGL). Part of the 
apparent noise in the signal is an effect of the sampling rate used to 
make the plot and a high frequency component in the sensor's output. The 
sensor is mounted on a nose-boom and a higher sampling rate clearly shows 
the boom's vibrational frequency in the signal. The radio altimeter is 
used only below 137m (450 ft) AGL; hence, only that portion of the data 
is shown in the figure. The radio altimeter has a known bias such that 
it reads 7.0m when the aircraft is sitting on the ground. 

Figure 5c shows the outputs from the triad of body-mounted accelero- 
meters and the outputs from the LTN-51 inertial navigation system (INS) . 

The body-mounted accelerometer signals are clearly very noisy, mostly due 
to the vibration of the aircraft's structure. The vertical accelerometer 
reading is compensated for the one-g reading it would normally have in 
steady, level flight before it is telemetered to the ground. The LTN-51 
output signals are much smoother because the unit has a "sum AV" mecha- 
nization; l.e., the instrument sums increments of velocity in discrete 
counts, and acceleration is inferred by differencing the sum at the be- 
ginning and end of a time interval and dividing by the length of the 
interval. For the plots the interval was 0,5 sec, so the result is effec- 
tively the average acceleration for each half second. The filters use an 
Interval of .05 seconds; nonetheless, the INS signals are much smoother 
than the sampled analog signals of the body-mounted accelerometers. Be- 
cause the INS uses a stabilized platform, the x and y acceleration outputs 
clearly show the two turns of the helix. 
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r Figure 8 presents the time histories of residuals for the MLS, TACAN, 

baromet *ic altitude and radio altitude measurements shovm on Fig. 5. These 
residuals are differences between the measurements and computed measurements 
based on the aircraft position derived from the ground tracking radar data. 
One cannot look at these residuals as true measurement errors since the 
radar tracking system has errors. Hence, if the plots do not behave as 
expected, the fault may lie in the tracking radar data as well as in the 
MLS, TACAN, or other measurements. 

The MLS range residual has an average value of about -12.2m (-40 ft) 
with high frequency noise of perhaps 3m rms. The bias and lower frequency 
component of the residual could be in the MLS measurement or it could be 
caused by radar tracking errors. 

The average MLS azimuth residual is -0.25 degrees. High frequency 
noise and very small anomalies occur in this residual at the same time as 
they occur in the actual recorded measurement (Fig. 5). Therefore, the 
noise and anamolies are indeed present in the MLS azimuth measurement and 
are not caused by any radar tracking errors. The bias, however, may pos- 
sibly be attributable to the radar tracking data. 

The MLS elevation data are obviously very bad until about 140 seconds. 
The residual at about 200 seconds is also due to an error in the MLS 
measurement. It is known that this residual is not due to an error in the 
tracking radar data because, if it were, a similar residual would occur 
in the barometric altimeter data. 

The TACAN range residual shows an apparent bias of about 107m (350 ft) 
with high frequency errors having an rms of about 30ra (100 ft). The TACAN 
bearing residual shows these data have little utility for nearly the com- 
plete helix. During times when these data are good, the scale is too 
large to conclude anything with regard to the residual characteristics. 

The barometric altitude residual shows characteristics of bias and scale 
factor errors. The bias starts at about -23m (-75 ft) at 805m (2640 ft) 


1 
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Figure 8. Residual Time Histories for On-board Measurements 










altitude MSL and reduces to around -9m (-30 ft) at 427m (1400 ft) altitude. 
The high frequency noise has an rms of about 6m (20 ft). The radio alti- 
tude residual indicates a scale-factor error and some anomalies before 
300 sec which are believed to be caused by the terrain. The spike In the 
residual at about 305 seconds can be seen In the barometric altimeter 
residual and the MLS elevation residual. Therefore, It Is assumed that 
this spike Is not caused by the radio altimeter but by bad radar tracking 
data at this time. The scale-factor error may be caused by any of the 
following: 

1. scale-factor error In the radio altimeter, 

2. slope of terrain the aircraft Is flying over, or 

3. radar tracking errors. 

It Is not known at this time which source Is the primary contrlbuter. 
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NAVIGATION SYSTEMS DESCRIPTION 


Overview 

All the navigation systems described here provide the estimated posi* 
tion and velocity of the VTOL aircraft by combining inertial measurements 
with measurements from the navaids. Figure 9 shows how the navigation 
systems are implemented for tests in the V/STOLAND avionics system. The 
basic computer contains all the primary software for operating the 
V/STOLAND system. The research computer is strictly for research 
investigations and the researcher supplies software for the particular 
experiment under test. The basic software is designed such that many different 
types of research experiments can be conducted without changing the basic program. 

As shown on Fig. 9 all data for navigation experiments except the 
LTN-51 INS accelerometer outputs come into the basic computer. All data in- 
put to the basic computer are also sent to the research computer. The 
switches shown on the figure are under the pilot's control. As can be seen, 
he may either use complementary filter or Kalman filter state estimates for 
driving the basic computer's display, guidance and control logic. The re- 
search-mode button controls vdiich state estimates are used. Also by use of 
keyboard inputs, he may select either the strapped down IMU or the LTN-51 as 
the source of acceleraclon input to the complementary filter in the basic 
computer or the Kalman filter in the research computer. 

Figure 10 is a block diagram illustrating the general structure and 
functions of all the navigation systems. The inertial measurement unit 
(IMU) provides sufficient data for calculating the aircraft acceleration 
in a runway referenced coordinate frame. The accelerations are integrated 
to keep the position and velocity estimates current. When hardware dis- 
cretes indicate the navald measurements are valid, their values are 
compared with estimated position data. If the difference satisfies the 
data rejection algorithm, then state corrections are calculated by a 
specified algorithm and added appropriately to the estimated state. 
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agram of Navigation Systems 


The vertical channel la handled Independently of the level channels in 
the systems described. The vertical channel is started using the raw 
barometric altitude reading for the vertical position at the initialization 
time point • and zero initial value is placed on the vertical velocity. 

When MLS elevation data become available it is used as the primary reference 
until the aircraft gets below about 152m. At this point the radio altimeter 
measurements become the primary vertical position reference. 

For the level channel, x-y position initialization is performed using 
MLS range and azimuth if available; otherwise, the less accurate TACAN range 
and bearing measurements are used. Airspeed and aircraft heading measure- 
ments are used to initialize the level components of velocity. 

The automatic measurements selection algorithm for the level channels 
will use MLS range and azimuth if available; otherwise, TACAN measurements 
are used. If neither source of data is available, a dead-reckoning mode 
involving either inertial information only or inertial information and 
airspeed measurements is used. 

The navaid reference changes as the aircraft enters the terminal area 
and precedes to the landing. During the transition from one reference to 
the next, transients in the estimated state occur. The block in Fig. 10 
called "navaid transition smoothing" is used in the Kalman filter to pre- 
vent these transients from causing rapid aircraft maneuvers. The transi- 
tion-smoothing logic does not exist in the complementary filters. 
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Cotcplenentary Filcers 


The complementary flltera which are used in the flight teata of the 
V/STOLAND avionics system were initially developed by Sperry Flight Systems 
(Ref. 2). The availability of some new data sources (MLS and INS accelero- 
meter data) a.id information gleaned from flight test results have led to 
a number of modifications. The complementary filter as currently 
mechanized in the Sperry 1819B basic computer is summarized in this section. 

Figure 11 is a block diagram of the complementary filter used in the 
V/STOLAND system. The MLS range, azimuth, and elevation, the TACAN range 
and bearing and the VOR/DME range and bearing measurements are fed through 
first-order pre-filters. In order to pre- "t lags caused by the pre-filter 
time constants, the estimated rates for each of the measurements based on 
the current state estimate are also fed to the pre-filter. Referenc-i selec- 
tion logic either manually through push buttons or automatically (if the 
auto nav mode is selected) determines which pre-filter navaid data are 
used for the raw x-y calculations. The rau x-y data and the acceleration 
in the runway reference frame as calculated from the raw inertial data are 
fed to the two third-order x-y navigation filters. The acceleration input 
source may be either the strapped-down IMiJ or the LTN-51, as was shown 
in Fig. 9. 

In the vertical channel the pre-filtered MLS elevation data and raw 
barometric and rad^o altimeter altitude data are fed to the reference 
selection logic. The reference selection logic is fully automatic for the 
vertical channel. The barometric altitude is used until MLS elevation data 
are valid. This is followed by a blending period where MLS elevation and 

estimated x-y coordinates are used to calculate one source of altitude and 
the barometric altitude is the second source. The raw altitude used for 
Che filter is a linear combination of the two altitude sources where the 
weight shifts with time from all barometric data to all MLS altitude data. 
The total time for the blend is 60 seconds. As radio altimeter data become 
available . radio altitude and the other source of altitude (biased baro. 
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Figure 11. *»IocV Diagr.i-; of Conpleaientary Filters 






blended baro-MLS, or MLS only) are also blended. In this instance the 
blending favors radio altitude as the altitude decreases. At 61m altitude 
the blending ceases and the radio altimeter Is the source. The raw 
altitude and the vertical acceleration are fed to a third order navigation 
filter for the altitude channel. 

Figure 12 is a block diagram of the pre-filter., osed in the comple- 
mentary filter. The filtered measurement is subtracted froc the raw 
measurement and the difference sent to a tolerance test. If the tolerance 
is exceeded, the raw measurement is rejected. If the tolerance test is 
passed, then the error signsl is limited before being multiplied by the 
reciprocal of the time constan: and integrated. The estimate' ’■.^te for 
the measurement is fed directl ' to the integrator for the filterec in 'asure- 
ment. The table on Fig. 12 gives the tolerance, limit level and time 
constants used for the pre-filters in the V/STOLAND complementary navigation 
filters . 

Figure 13 shows the third-order navigation filter for the x channel of 
the complementary filter. The y channel is identical in structure and filter 
gains. The switches in the figure are shown in the normal "navigation valid" 
operation of the system. In this mode of operation the estimated position 
x^ is subtracted from the raw position and the difference used as feed- 
back through gains and into the three integrators of the filter. 

The measured acceleration from the selected source feeds the integrator 
whose output is the estimated velocity Xj^. In this mode of operation the 
values of the filter gains depend on the source of the navaid-derived position. 
It should be noted that the pre-filter time constants shown of Fig. 12 also 
depend on the source of the navaid-derived position. The combination of 
the filter gains and the pre-filter time constants were selected so that, 
when MLS is in use, the overall complementary filter is more responsive in 
tracking the navaid-derived position than when TACAN or VOR/DME is in use. 

It should be mentioned that the pre-filters used in the V/STOLAND navigation 
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Figure 12 . Block Diagram and Parameters of the Prefilters 
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Figure 13. Third-Order Navigation Filter for x Channel 












system cause coupling to exist among all the channels. As a result • the 
stability of the navigation system is difficult to analyze, and simulated 
results were used to select the gains. 

In the normal mode of operation the components of wind are estimated 
in a runway referenced coordinate frame. This is achieved by sending the 
difference between measured airspeed and ground speed Xj^ into a first 
order filter with a 100 second time constant. In the dead reckoning mode 
the wind estimate is frozen and used with the airspeed and ground speed 
data to stabllze the ground velocity. In this instance the estimated 
ground velocity is stabilized at the airspeed value compensated with the 
last value of estimated wind. The wind estimate is not used for any other 
function. 

Figure 1^ shows the third-order navigation filter for the altitude 
channel. In this instance the configuration and gains are not dependent 
on the source of altitude data. When either the barometric or radio alti- 
meter is the altitude source, the filter dynamics are as given in the 
figure. As a result of the pre-filtering of MLS data and the blending 
algorithms, the filter's dynamics change substantially during an approach. 

Comments on the V/STOLAND Complementary Filter 

The V/STOLAND complementary filters navigation system has some un- 
desirable characteristics which should be removed in an operational design. 

1. There is no provision for providing smooth transition from one 
navaid source to another source. If the aircraft is in the 
automatic flight mode, undesirable steering transients occur when 
navigation aids are changed (for example, a transient occurs 
when switching from TACAN to MLS) . 

2. The pre-filters used in the V/STOLAND system introduce much 
unnecessary complexity without providing any improvements 

in performance. There is some evidence from simulation studies 
that Indicates the performance may actually be degraded. 
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Figure 14. Third-Order Navigation Filter for Altitude Channel 






3. The filter is too responsive to MLS measurement errors when 

the aircraft reaches hover. This leads to excessive undesirable 
maneuvering during the letdown phase. 

4. The navigation errors are relatively large during the helix 
descent and a moderate time thereafter. The filter is not 
responsive enough in tracking the MLS data in these regions. 


Kalman Filter 

This section describes the Kalman filter navigation system developed 
for the V/STOLAND avionics system. The navigation system is a modified 
version of the airborne software originally developed and tested in the 
STOLAND system (Ref. 3). 

Figure 15 is a block diagram of the V/STOLAND Kalman filter navigation 
system which is implemented in the research computer. The raw inertial 
measurement unit (IMU) data feed the block labeled "acceleration calculations". 
The calculations performed in this block depend on which source (strapped 
down or INS) of inertial data is used and the selected reference frame for 
acceleration bias estimates. The runway-referenced accelerations obtained 
from the acceleration calculations drive the navigation equations, which in 
turn provide the estimated positions and velocities. These two functions 
operate at 20 Hz in the mechanization. 

Raw navaid data feed the block labeled "measurement rejection and 
preprocessing". Here, hardware vallds and software tolerances are used to 
establish validity of the measurements. For valid measurements, the differences 
between raw measurements and measurements computed from the estimated state 
(called residuals) and their associated partlals are formed and accumulated. This 
logic is executed at 10 Hz. Once per second, measurement selection logic 
picks the desired set of residuals and partlals and sends them to the filter 
algorithms. The filter algorithm executes a square-root implementation of 
the Kalman filter and produces the estimated error state. The estimated 
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error state is used to correct the position and velocity estimates, the 
acceleration bias estimates, and the TACAN range and bearing bias estimates. 
It also feeds the smoothing logic along with the estimated state. The 
smoothing logic prevents abrupt changes in the position and velocity which 
are sent to the basic computer for display, guidance and control purposes. 

The filter algorithm is actually two Independent implementations consisting 
of an x-y filter and a z filter. The level channels (x-y filter) use the 
MLS range and azimuth or the TACAN range and bearing to estimate the 8-ele- 
ment error states given by 

dx(l) * error in x component of position, (dx) 

dx(2) • error in y component of position, (dy) 

• 

dx(3) “ error in x component of velocity, (dx) 

dx(4) “ error in y component of velocity, (dy) 

dx(5) * error in 1st component of level acceleration bias, (b^^) 

dx(6) * error in 2nd component of level acceleration bias, (b ) 

ay 

dx(7) ■ error in TACAN range bias estimate, (b^) 
dx(8) * error in TACAN bearing bias estimate, (b^) 

The acceleration bias estimates can be in either the runway reference frame 
or a path-referenced frame depending on pilot input via the keyboard. 

The implementation has both automatic measurement selection and manual (by 
means of push buttons) modes. In the automatic mode, MLS will always be 
used if both range and azimuth residuals are available from the preprocessing 
logic. 

The vertical channel (z-f liter) uses the barometric altitude, altitude 
computed from MLS elevation or radio altimeter altitude to estimate the 
3-element error state given by, 

dz(l) ■ error in z component of position, (dz) 
dz(2) ■ error in z component of velocity, (dz) 
dz(3) ■ error in vertical acceleration bias, (b^^) 
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For Che z channel the measurement selection algorithm Is automatic only. 
Radio altitude Is used when the radio altimeter output Indicates the alti- 
tude Is less than 137m and the hardware flag Indicates valid data. MLS 
elevation data will be used otherwise If hardware and software tests show 
it is valid. Barometric altitude Is used In all other cases. A special 
algorithm is used to estlti>ate the bias In barometric altitude and to decay 
this estimated bias, as will be explained. 

The overall Kalman filter mechanization for the V/STOLAND system Is 
very similar to the STOLAND system dlscrlbed in Ref. 3. The next two 
sections describe the areas where differences exist. The reader should 
consult Ref. 3 for those areas not covered here. 

Level (x-y) Filter 

The V/STOLAND x-y filter differs from the STOLAND x-y filter in the 
following areas. 

1. MLS range and azimuth are used instead of the MODILS range and 
azimuth. 

2. Airspeed Is only used In Che initialization of the x and y 
components of velocity. Winds are not estimated in this 
system. 

3. The estimated acceleration biases can either be In the runway 
frame or in a path-referenced frame. The option Is controlled 
by the pilot. 

4. The source for runway-referenced acceleration can be either 
the strapped-down IMI3 or the LTN-51 INS. This option is also 
selected through the keyboard. 

5. Initialization mode uses the MLS data, if available; otherwise, 
TACAN data are used. 
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The transition matrix, 4>, is approximated as 

« I + A (1) 

The non-zero elements of A in £q. (1) for the x-y channel are given by 


A^(l,3) 

A^(2,4) - At 

A^(l,5) 

A^(2,6) -(cos i>)ht^/2 

A^(l,6) 

-(sin i^)At^/2 


(sin 

A,(3,5) 

A^(4,6) -(cos \J^)At 

A^(3.6) 

-(sin4;)At 

A^(4.3) 

(sinti^)At 

A^(5.5) 

A^(6,6) —At/t^ 

A^(7,7) 

-At/x 

r 

A^(8.8) 



Here, 

At ■ period over which transition matrix is used 

T * • time constant for acceleration colored noise (20 sec) 

d 

» time constant for TACAN range colored noise (1000 sec) 

T^* ■ time constant for TACAN bearing colored noise (1000 sec) 

When the path referenced system is used, the of Eq. (2) is the air- 
craft heading with respect to the runway. In the other options "'ii of Eq. 

(2) is set zero. The path-referenced system is only used with the strapped- 
down INO and was developed to determine if this implementation would alleviate 
navigation offsets following helix exit. 


* These variables are input to the program as reciprocals and may be 
modified by pilot Inputs through the keyboard in even Increments of 
IZ from 0 - 400%. The value shown is the 100% value. 
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The non-zero elements of the forcing matrix 4>u^ (see Eq. (A. 4) of 
Ref. (3)) are given as 


^^(3,3) 

♦u^(5,5) 

4>u^(7,7) 

^^(8,8) 


/IT 


^u^(6,6) - 0^ sTT^TT 


°r \/ 2At/ij. 

- v^' 2it/x 






( 3 ) 


here 


At 

0 

V 

a 

a 

0 


r 



period of the covariance matrix update (1 sec) 

standard deviation (std) of velocity noise (.0762 m/s) 

2 

std of acceleration colored noise (.1524 m/s )* 
std of TACAN range colored noise (304.3 m) * 
std of TACAN bearing colored noise (2 deg) * 


The runway referenced accelerations are given by 
X “ XriLni b COSlp ~ b siWp 

IMU ax ay 


where ip In Eq. (4) is the aircraft heading with respect to the runway 
when the path-referenced system is used and zero otherwise. 


MLS range and azimuth - The MLS measurements used in the x-y portion 
of the Kalman filter are range and azimuth from a co-located DME transponder 
and azimuth scanner. The range measurement is modeled as 


Here, 


Y 

mr 




(x-x )^ + (y-y + q 

m H} m tnr 


(5) 


X ,y , 
m "^m 


z 

m 




coordinates of the MLS transponder and scanner 
with respect to the runway reference frame, 
the random noise error in the range measurement. 


* The standard deviations may be modified by pilot Inputs through the 
keyboard in even Increments of 1% from 0-400X. The value shown Is 
the 100% value. 
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The estimated measurement is computed from 


Y - v/(x-x ) + (y-y_) + 

IDr U HI HI 


Here, x and y are state variables obtained from the x-y filter, and z 
is obtained from the z-filter. 


The non-zero elements of the row vector H for the range measure- 
ment are calculated from 




The variance of the random noise error in the range measurement is 
assumed to be a constant given by* 


0^^ - (18.3 m)' 


The MLS azimuth measurement is modeled as 


iDd n in iH& 


Here, is a random error in the azimuth measurement. 


The estimated measurement is computed from 


tan I-(y-y )/-(x-x )]. 

m m 


Again, x and y are state variables of the x-y filter, and z is obtained from 
Che z filter. 


* The standard deviations of the MLS range and azimuth are keyboard input 
quantities with range 0-400%. The values shown are 100% values. 
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the azimuch measure- 

( 11 ) 


x-y Filter Initialization - As mentioned previously, the V/STOLAND 
x-y filter will be initialized from MLS range and azimuth if both are valid 
at the time the initialization is requested. Otherwise, TACAN range and 
bearing are used if they are both valid. Initialization will not take place 
until one or the other source of raw x-y data is available. The TACAK 
initialization equations are given in Ref. 3. The MLS initialization 
equations for the state and the square-root covariance are given here. 


The non-zero elements of the row vector H for 
ment are given by 

max 'm m m 

H (2) - -(x-x )/t(x-x + (y-y )^1 

max n m m 


■ (.1 deg)' 


Here 


X ■ 

*m - 

r cos 
c. 

(Y ) 
ma 

y ■ 

^m 

r sin 
c 

ma 


(13) 



h “ altitude of aircraft above the MLS station computed 
ma 

from barometric altitude. 


The velocity components are computed from airspeed using equations given 
in Ref. 3. 


The non-zero elements of the Initial square root covariance matrix 
for the x-y filter are given for MLS initialization by the following: 
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1 


Wx(l.l) 

■ 

0 cos(Y -) 
mr ma' 

Q 

B 

■ 

W^(1.2) 

■ 

sln(Y„) 
mr ma 


Wx(2.1) 

■ 

0 (y -y) 
ma'^m ^ 

0 ■ 

ma 

W^(2,2) 

• 

„ (x -x) 
Oma m 


Wx(5,3) 

■ 

-cos(i^.-\|; , 

i r va 

0 ■ 

va 

Wj^(5,4) 

■ 

-sln(4'.-'/^ )‘^ , 
1 r va 


Wx(6.3) 

■ 


Q 

1 

Wjj(6,4) 

• 



W^(7,3) 

m 

0 t 

wx 

0 - 

wx 

W^(8,4) 

m 

0 

wy 

0 • 

wy 

W^(9,5) 


°ax 

m 

W^(10,6) 

m 

°ay 

U 

ax.v 


standard deviation (std) of 
MLS range noise (18. 3n) 

std of MLS azimuth noise (.1 deg) 


std of random error in air- 
speed measurement ~ 0.61 m/sec. 

std of bias error in initial 
h<sading ~ 2 deg. 

std of X component of wind 
^ 6.1 m/sec. 

std of y component of wind 
» 6<1 m/sec. 

acceleration random error 
std ■ . 3 ra/sec- 


2 Filter 

The 2 filter for the V/STOLAND system is a 3-state Kalman filter 
as was mentioned earlier. The barometric altitude bias is estimated in an 
ad-hoc manner external to the filter. The filter is Inltialixed from baro- 
metric altitude information; and the estimated bias in barometric altitude, 
estimated vertical acceleration bias and the estimated vertical velocity 
are set 2 ero. Following initialiration the operation of the filter is as 
described in the following. 


a) Only one source of altitude measurement is processed by the 

Kalman filter at a time. If (i) the radar altimeter valid is 
set, (ii) the estimated altitude above ground is less than 183m, 
and (ill) the radio altimeter altitude measurement is less than 
137m, then the radio altimeter measurement is processed. If (i) 
radio altimeter measurements fail the above tests, (11) the MLS 
elevation valid is set, (lii) MLS elevation is less than 10°, 
and (iv) MLS range and aximuth are being used by the level KLalman 
filter, then the MLS elevation measurement is processed. If 
both radio altitude and MLS elevation are not available, then the 
barometric altitude measurement adjusted for the current baro 
altitude bias, if any, is processed. 
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b) The barometric altimeter bias Is estimated apart from the Kalman 
filter logic. The Initial estimate of the bias Is zero. In 
the first Kalman cycle wherein MLS elevation data or radio alti- 
meter data are accepted according to the criteria stated In 
(a) above, the baro bias is estimated. The accepted altitude 
measurements (radio altimeter or MLS elevation) are averaged 
and any measurement differing from the average by more than 6.1m 
is discarded. If at least 6 measurements remain, the average is 
recomputed, and the baro bias is estimated as the difference 
between this average and the current baro altimeter measurement. 


c) 


This estimate is used to initialize a first-order filter which 


operates as long as the MLS elevation or radio altimeter measure- 


ments are accepted. This filter is 

hbb<t) 

+ d-e”*’^) (h + z) 
m 


(14) 


Here 

/N 

h,, is the baro bias estimate, 
bb 

h is the altitude above ground measured by 

the radio altimeter or MLS elevation (positive up), 

z is the aircraft altitude estimate in the Kalman 
filter (positive down, referenced to runway). 


At is .05 sec. 


d) If both the radar altimeter and MLS elevation data are rejected 
according to the criteria of (a) above, then the baro altitude 
bias estimate is held constant at the last value computed by 
the filter in (c) above as long as the aircraft descends. If 
the aircraft ascends, the bias estimate is reduced to zero. 
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This is done in ateps of l/8th of its last filtered value for 
each 61in gained above the altitude where the filter ceased, 
until the absolute value of the remaining bias is less than 2,4m. 
This 2,4m tolerance avoids possible difficulties due to com- 
puter truncation errors and is small enough that it will not 
affect navigation performance. For example, suppose the bias 
filter has a current estimate of 16m when the MLS elevation 
data which were driving the filter become unacceptable, and 
suppose this occurs at 300m altitude. As long as the air- 
craft descends, the bias remains at 16m. If the aircraft 
reached 180m and then climbs, the bare bias estimate will be 
reduced to 14m at 361m altitude, 12m at 422m altitude, etc. 

At 727m the bias estimate is reduced to 2ra, and no further 
reduction will take place. 

e) If the MLS or radar data again become acceptable within 20 
seconds, the bias filter resumes operation at the then pre- 
sent bias estimate. If more than 20 seconds elapses, the 
initialization procedure described in (b) above is repeated 
and the filter begins operation with that bias value. 

f) If the bias filter is being driven by MLS elevation data when 
radio altimeter data become acceptable, then (i) the bias 
filter is reinitialized according to the procedure in (b) 
using the radio altimeter data, (11) the bias filter is 
driven by the radio data alone, and (Hi) the MLS elevation 
data, although acceptable, are ignored. 

The 3-state z filter and ad-hoc baro bias algorithm were developed 
for V/STOLAND for the following reasons. 

(1) The aircraft should fly baro-referenced altitude in the terminal 
area, even if a good source of true altitude is known, because 
it is assumed other traffic is also using baro altitude. 

(2) As the aircraft descends to land, true altitude must be the 
desired goal of the navigation system. Should MLS fail or not 

be available during part of the approach, then the best procedure 
for estimating altitude is the baro altitude measurement adjusted 
for the last calibration of the bias. 
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(3) In case of go-around the system reference must go back to the 
unbiased baro reference, since this remains the reference used 
by other aircraft. 

(4) The barometric bias estimates and z filter performance were not 
acceptable when using the 4-state Kalman filter of Ref. 3. The 
problem was traced to modeling inaccuracies combined with the 
fact that the baro bias state is unobservable except when another 
source of altitude is available. These factors led to very poor 
estimates. 

Note that the actual error in the barometric altimeter is more accurately 
modeled as a scale factor error and a bias error (see Fig. 8). Further 
effort could perhaps lead to a better model for the barometric altitude error 
and performance improvements over the 3-state filter described here. 

MLS Elevation - Define the aircraft relative position coordinates 
with respect to the MLS elevation antenna 



(15) 


Here (x„,y„,x„) give the position of the MLS elevation antenna in the runway 
c c r- 

reference frame. 


Let 



(16) 


r 

V X, 




(17) 


The altitude measurement calculated from the MLS elevation measuremint is 
expressed as 

Y “ r,tan(E:) + q 
e 1 ' ^e 
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Here 


c ■ the MLS elevation measurement • 

■ the random noise error in the pseudo-altitude measurement . 

The estimated measurement is given by 

Y - -z (19) 

e 

The non-zero element of the row vector H is given by 

H (1) = -1 (20) 

ez 

The variance of the random noise error in the measurement is assumed to be 

the range-dependent quantity 

Q = (0 r )^ (21) 

e e e 

Here 

= std of MLS elevation noise (.1 deg*) > 
r^ = calculated ranged to MLS antenna, Eq. (16). 


Keyboard Controlled Quantities 

The V/STOLAND system is designed with a pilot-operated keyboard for 
changing program quantities. Table 1 shows the Kalman filter quantities 
which may be changed by keyboard entry during flight tests. The software 
is arranged so the pilot enters the desired percentage of a nominal value 
The software checks to see if this percentage lies within the range shown 
in the table and, if so, modifies the quantity to the desired percentage 
If the input lies outside the allowed range, then the input is limited to 
lie within the range, and the change is given to the quantity. 


* The std of the elevation may be modified by keyboard input. 


54 


Table 1 


Quantities Modifiable through the Keyboard 


Quantity 

Symbol 

Nominal (100%) Value 

Range (%) 

std acceleration 
noise x channel 

a 

ax 

.15m/s^ 

O-AOO 

std acceleration 
noise y channel 

0 

ay 

.15m/s^ 

0-400 

reciprocal of 
time constant 
X channel 

1/t 

ax 

.05sec ^ 

0-400 

reciprocal of 
time constant 
y channel 

1/t 

ay 

.05sec ^ 

0-400 

Initial std of 
acceleration error 
x-y channel 


. 15m/s^ 

0-400 

MLS range 
noise std 

0 

mr 

18,3m 

0-400 

MLS azimuth 
noise std 

o 

ma 

• 2deg 

0-400 

TACAN range 
noise std 

Q 

ft 

183m 

0-400 

TACAN bearing 
noise std 

°tb 

2deg 

0-400 

Barometric alti- 
tude noise std 

°hb 

60m 

0-400 

MLS elevation 
noise std 

0 

me 

• 2deg 

0-400 

Radio altimeter 
noise std 

0 

r 1 

. 3m 

0-400 

std acceleration 
noise z channel 

0 

az 

. 15m/ s^ 

0-400 

reciprocal of 
time constant 
z channel 

1/t 

az 

.OOlsec ^ 

0-400 

Option flag for 
acceleration 


100 

+100 = strapped-down 
-100 = LTN-51 INS 

source 




Option flag for 
acceleration bias 
frame 

“ 

100 

+100 = runway reference 
-100 >= path reference 
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Smoothing Logic 


The smoothing logic used In the V/STOLAND system Is very close to that 
used In STOLAND (see page 114-115 Ref. 3). The only significant difference 
Is In the decay algorithm for the position and velocity smoothing vector. 

In Ref. 3 . this algorithm Is given as 


c (t + At.) ■ e (t) 

X f X 

c (t + At,) = (t) 

VI V 


( 22 ) 


In the V/STOLAND system the change In c^ and c^ given by Eq . (23) 


Is formed as 


Ac^(t + A t^) = (1 - e~-°^^'^^)c^(t) 
Ac^(t + A t^) - (1 - e'*°^^^v)c^(t) 


(23) 


If the absolute value of the change Is less than a desired Input value , 

then Eq. (22) Is used. If the change is larger than the limit value, then 

the limit value is used with the appropriate sign to decrease the magnitude 

of c and c . 

X v 

Current limit values used in the V/ STOLAND system are: 

L ■ 3m/s for x-y position changes, 

^ » y 

* 1.5m/s for z position change, 

L« • • “ .61m/s/s for velocity changes. 

X , y , z 

The time constants T and are currently set at 10 seconds. 

XV 

This algorithm serves the purpose of smoothing the once-per-second 
changes calculated by the Kalman filter algorithm as well as smoothing the 
transition when navald sources are switched. For example, the large position 
change of perhaps 100m which occurs when switching from TACAN to MLS will 
require about 33 seconds to get into the position estimate which drives 
the guidance display and control logic. The change is slow enough that 
only a very minimal, smoothed transient in aircraft maneuvering occurs. 
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CONCLUDING REMARKS 


Navigation systems for research investigation in the V/STOLAND avionics 
system have been described, These systems are currently being studied in 
flight tests to evaluate their performance, and further refinements are 
likely. 

The general problem of operating a VTOL aircraft from a heliport in 
all-weather conditions appears to require higher accuracy and smoothness 
of the estimated state than for a CTOL aircraft. These requirements arise 
in part due to the small size of the heliport. In addition, lateral 
maneuvers of the aircraft during the final letdown phase are intolerable 
to most pilots. 

Flight tests of these systems will provide information on whether a 
single MLS system at an airport terminal can satisfy both VTOL and CTOL 
aircraft landing navigation requirements as well a^ information on the 
requirements of the onboard inertial sensing and software sophistication 
for the VTOL system. 
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